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Abstract 

Utilizing the unique and rehable uhrasmah-x predictions of the dynamical (ra- 
diative) parton model, nominal event rates and their detailed energy dependence 
caused by a variety of cosmic UHE neutrino fluxes are calculated and analyzed. In 
addition, maximal Regge-model inspired small-x structure functions are employed 
for obtaining optimal rates which do not necessarily require 'new' physics interpre- 
tations. Upward //^ + /x~ event rates are estimated by taking into account total and 
nadir-angle dependent regeneration effects due to neutral current interactions. For 
exploring extragalactic neutrino sources at highest energies 10^ GeV) with mod- 
ern (future) ground-level telescopes, we analyze horizontal air shower event rates 
and shower events caused by Earth-skimming tau-neutrinos, in particular their de- 
tailed shower- and cosmic neutrino-energy dependence. As an illustration of 'new 
physics' implications we estimate the relevant horizontal air shower event rates due 
to spin-2 Kaluza-Klein 'graviton' exchanges in neutral current neutrino-quark and 
neutrino-gluon interactions at low TeV-scales. 



1 Introduction 



Detection of cosmic ultrahigh-energy (UHE) neutrinos with energies above 10^^ eV is one 
of the important challenges of cosmic ray detectors in order to probe the faintest regions 
of the Universe that are otherwise shielded from us by large amounts of matter. Their 
observation will probe particle (possibly 'new') physics as well as astrophysics phenomena 
such as galaxy formation. The sources of UHE neutrinos range, however, from the well 
established to the highly speculative [1, 2, 3, 4, 5]. The well measured cosmic rays up to 
the GZK cutoff [6] at around 5 x 10^^ eV, where they necessarily interact with the 2.7 K 
cosmic microwave background through p7 mr'^, produce the 'guaranteed' (cosmogenic) 
flux of 'Greisen neutrinos' when the pions decay [6, 7]. In addition, far larger neutrino 
fluxes are predicted in models of active galactic nuclei (AGN) [8, 9, 10], gamma ray bursts 
(GRB) [11], decays of exotic heavy particles of generic top-down or topological defects 
(TD) [12, 13, 14, 15, 16] and Z-bursts [17, 18]. Representative fluxes of some hypothesized 
sources are displayed in Fig. 1 which we shall use for all our subsequent calculations. For 
illustration the steeply falling background atmospheric (ATM) neutrino flux [19, 20, 21] is 
shown as well which originates from (anti)neutrinos produced by cosmic ray interactions 
in the Earth's atmosphere. 

Besides these violently different model expectations for ( ant i) neutrino fluxes, there are 
further uncertainties when calculating event rates for neutrino telescopes due to the sen- 
sitivity of ^I'^N cross sections to the parton distributions in the yet unmeasured ultrasmall 
Bjorken-x region, x < 10~^, not accessible at the weak scale Mw by deep inelastic scatter- 
ing (DIS) experiments. At the highest (anti) neutrino energies E,^ shown in Fig. 1, contri- 
butions from the region around x ~ M^/2M]^Ei, ~ 10~^ to 10~^ become non-negligible 
and therefore known structure functions and parton distributions at x ^ 10~^ have to be 
extrapolated to a; < 10~^ as soon as E,^ 10^ GcV. (Here one commonly assumes that the 
DIS small-x measurements at x ^ 10~^ and moderate momentum scales can be safely 
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evoluted to the relevant scale — by standard QCD renormalization group (RG) 
techniques.) Such extensive extrapolations are performed either of fits to existing data 
using specific, possibly arbitrary and unreliable, assumptions (e.g. various 
fixed power behaviors in x of structure functions as a; — 0) [20, 21, 22, 23, 24, 25, 26] or 
by using [27, 28] the QCD inspired dynamical (radiative) parton model which proved to 
provide rehable high energy predictions [29, 30, 31, 32] in the past [33, 34]. Within this 
latter approach the entire partonic (gluons and (anti)quarks) structure at small fractional 
momenta x <, 10"^ can be understood and calculated via RG evolutions from first princi- 
ples, i.e. QCD dynamics, independently of any free (fit) parameter in the small-x region 
due to valence-like gluon and sea input densities at some low momentum scale Qo ~ 0.5 
-0.6 GeV [30, 31, 32]. Having successfully predicted the small-x behavior of structure 
functions between x — 10~^ to 10"^ prior to experiments, it is not unreasonable to expect 
the unambiguous dynamical results of the radiative parton model between x — 10~^ and 
10^*^ or 10~^ to be reliable as well. These ultrasmall-a; predictions are furthermore per- 
turbatively stable and unique at the relevant momentum scales ~ [32] . Moreover, 
consistent BFKL model resummations of subleading In ^ contributions in leading order 
QCD yield remarkably similar results [35] even at highest neutrino energies of 10^^ GeV. 

A detailed analysis [36] has shown that within the radiative parton model all relevant 
(anti)neutrino-nucleon cross sections can be calculated with an uncertainty of typically 
about ±20% at highest neutrino energies of 10^^ GeV. We shall adopt this approach for all 
our subsequent calculations, in particular the canonical GRV98 parton distributions [32] 
in next-to-leading order QCD will serve as our appropriate nominal benchmark set of den- 
sities. Notice that cross sections obtained from the fitted CTEQ3-DIS parametrizations 
[37] Bit X ^ 10~^ with their assumed fixed-power extrapolation to x < 10~^, as used in [20], 
accidentally coincide practically with the ones derived from the ultrasmall-x dynamical 
predictions of the radiative parton model [36]; only at highest neutrino energies of 10^^ 
GeV, i.e. x ~ 10~^ to 10~^, they are about 10% larger than the radiatively generated ones. 
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Actual calculations of event rates could therefore also be performed with the CTEQ3- 
DIS densities which are easier to use since heavy quarks (c, b, t) are effectively (but less 
adequately) treated as massless intrinsic partons, in contrast to the more appropriate but 
also more cumbersome explicit calculation of massive heavy quark contributions required 
when using the 'fixed flavor' GRV98 densities [32]. (A comparison of the two approaches 
can be found in [36].) Similar remarks hold true for the CTEQ4-DIS parametrizations 
[38], used in [21], which underestimate, relatively to GRV98, the neutrino cross sections 
by about 20% at ~ 10^° - lO^^ Q^y 

The detection of UHE cosmic neutrinos is, however, extremely difficult. For ener- 
gies below 10® GeV it is widely believed that one of the most appropriate techniques for 
neutrino detection consists of detecting the Cerenkov light from muons or showers pro- 
duced by CC and NC neutrino interactions of mainly upward-going neutrinos in large- 
volume underground water or ice detectors such as AMANDA/IceCube, NESTOR and 
ANTARES [2, 39]. However, above 40 TeV the Earth's diameter exceeds the interac- 
tion length of neutrinos and their shadowing in the Earth rapidly increases above 100 
TeV [20, 21] which severely restricts rates in underground detectors which are bounded 
by detection volumes of at most 1 km^. Eventually it becomes beneficial to look for 
events induced by downward-going and (quasi) horizontal neutrinos, provided of course 
that downward-going events produced by interactions within the instrumented under- 
ground detector volume can be efficiently observed. At energies above 10^ GeV, where 
the (anti) neutrino interaction length is below about 10^ km water equivalent in rock, 
upward-going neutrinos are blocked by Earth and thus under-water and under-ice km- 
scale detectors become ineffective due to the opaqueness of Earth to upward neutrinos. 
Therefore large-area ground arrays or surface fiuorescence telescopes such as AG AS A, the 
HiRes detector, the Pierre Auger Observatory and the Telescope Array [26, 28, 40, 41], 
for which the interaction medium is not the Earth but the atmosphere, could become 
instrumental in exploring the whole spectrum of cosmic neutrino fiuxes up to highest 
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neutrino energies of about 10^^ eV shown in Fig. 1. In particular, a novel alternative 
method for detecting UHE ( ant i) neutrinos with > 10^ GeV has been recently sug- 
gested via Earth-skimming neutrinos [42, 43, 25], dominantly z/^ travelling at large nadir 
angles producing 'double bang' events in the atmosphere when converting to tau-leptons 
exiting the Earth's surface which, when decaying, produce the second shower bang. 

When upward-going UHE neutrinos penetrate through the Earth, they undergo atten- 
uation (absorption) due to charged and neutral current interactions as well as regeneration 
[44, 45] due to the neutral current interactions (which shift their energy rather than ab- 
sorbing them) at high energies. The NC shifts the neutrino energies to lower energies 
and, in addition, regeneration populates the lower energy part of the flux spectra shown 
in Fig. 1. In Sect. 2 we shall perform detailed calculations of these regeneration effects 
[35] not only for our representative fluxes in Fig. 1 but in particular also for the expected 
event rates for modern underground detectors. Depending on the shape of the neutrino 
flux, these effects can be substantial, i.e. increase the non-regenerated event rates on the 
average by about 20%. For comparison we therefore have also to recalculate, for our rep- 
resentative fluxes in Fig. 1, the upgoing event rates where regeneration effects have not 
been taken into account when calculating the shadow factor [20, 21, 23]. In addition we 
show in detail what high-energy bins of the original fluxes in Fig. 1 remain experimen- 
tally accessible after their depletion and regeneration due to NC interactions at highest 
neutrino energies. Furthermore, downward events initiated by UHE neutrinos will be 
also calculated using recent Regge-inspired power-like small-x extrapolations in order to 
learn about possible upper bounds of event rates which may be accommodated by 'con- 
ventional' standard model approaches and which do not necessarily require 'new' physics 
interpretations. Along similar lines we study in Sect. 3 event rates of quasi-horizontal air 
showers to be detected by large area ground arrays and surface air fluorescence telescopes, 
as well as of Earth-skimming UHE i/^-'s which will allow to test cosmic neutrino flux mod- 
els at highest energies, Ei, >, 10® GeV. We shall in particular be interested to what extent 
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shower rates as measured in certain bins of £^sh or Ej- will be able to explore and delineate 
specific features of the energy profile of the various predicted incoming cosmic neutrino 
fluxes ^{E,y) in Fig. 1. Finally, as an illustration of possible implications of 'new' physics, 
we present in Sect. 4 highly speculative expectations of string theories with large 'extra 
dimensions' for the relevant quasi-horizontal air shower event rates. Section 5 summarizes 
our results. 



2 Event Rates for Underground Cerenkov— Detectors 



The upward-muon event rate depends on the ^jy^N cross section through the interaction 
length Aint that governs the attenuation of the original neutrino flux due to interactions in 
the Earth as described by the shadow factor S as well as through the probabihty that 
the neutrino converts to a muon energetic enough to arrive at the detector with energy 
larger than the assumed threshold energy -B™"^. The total (nadir angle integrated) 
upward-muon event rate per second in a detector with an energy dependent area A{Efj^) 
is then given by [20, 21] 

-j^ ■gmin ^ /f) 

rate = 27r / dE, / ' ' dy A(E^)P^(E„ y; E--)S(E,) (2.1) 

where the factor of 27i is the effective solid angle for upward muons, i?^ = (1 — y)Ei, 
and ^,^{E,^) refers to the original cosmic neutrino fluxes shown in Fig. 1 incident on the 
surface of the Earth (i.e. = $/4) before they undergo attenuation and regeneration. 
Unless otherwise stated, the energy integration is always extended to the highest energy 
of 10^^ GeV in Fig. 1. Furthermore, 

P,(E^, y; Ef-) = R(E„ Ef^) (2.2) 

with A^^ = 6.022 x 10^^ and where all relevant expressions for calculating CC and NC 
^i/W cross sections in the 'fixed flavor' and 'variable flavor' scheme can be found in [36]. 
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It may be helpful to notice that, for an energy-independent effective detector area Aes, 
(2.1) reduces to the more common expression [20, 21] 

rate = 27rAeff / dE, P^(E,; E^-) S(E,) — ^ (2.3) 

with 

P,(E,; Ef^) = Na a'^^ciE.) {R {E.; Ef^^)) (2.4) 
and the average range {R) of a muon in rock is given by [20] 

{R{K;Ef^))^-J-- / dyR{{l-y)E.,E--) ^^cciyy) _ ^3.5) 

The range R of an energetic muon in (2.2) and (2.5) follows from the energy-loss relation 
[46] 

-dE^/dX = a^{E^) + (3^{E^) E^ (2.6) 

with X being the thickness of matter traversed by the muon in units of g/cw? = cm 
we. Despite the very weak energy dependence [46] of the ionization loss a^^En) we use 
[20, 46] = 2.0 x 10~^ GeV (cm we)~^ since the effect of a,j, is neghgible for high energies. 
Assuming the fractional energy loss /3^(£'^) to be energy independent as well, for example 
[20, 21, 46] (3^ = 3.9 x 10"^ (cm we)"^ or [25, 46] (5^ = 6.0 x 10"^ (cm we)-^ the relation 
(2.6) can be integrated analytically, 

R(E., E^'"') = X(E^''') - X(E.) = — In ^ i^C^^ • (2.7) 
Alternatively, using a QED-oriented energy dependence [47] 

a^{E^) = [2.033 + 0.077 ln(£;^/GeV)] x 10"^ GeV (cm we)~^ 

/5^(E^) = [2.229 + 0.2 ln(E^/GeV)] x 10"^ (cm we)"^ , (2.8) 

(2.6) has to be integrated numerically. In Fig. 2 we show the average muon ranges 
resulting from these three different choices and compare them with the Monte Carlo result 
of Lipari and Stanev [48] where regeneration has been taken into account and which has 
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been presented only up to 10^ GeV. Since this latter Monte Carlo result agrees best with 
our result using a constant and = 6.0 x 10~^ (cm we)~^, this choice appears to 
be best suited for extrapolations beyond 10^ GcV and we therefore will utilize it for our 
subsequent calculations. Needless to say that this conclusion is independent of the specific 
choice of parton distributions since obviously {R) in (2.5) does practically not depend on 
them. (If instead the numerical range of Fig. 2 were used, the resulting event rates would 
be about 30% larger.) 

Next we turn to the shadow factor in (2.1) describing the attenuation (due to (Jtot = 
ccc + ctnc) and regeneration (due to aNc) of ( ant i) neutrinos when penetrating through 
the Earth. Both effects are summarized in the transport equation [44, 45] for the neutrino 
fiux ^^{E,X) at 'depth' X 

= 1 + /" ^ ^<i^^^^E„X) (2.9) 

dA Aira,{J^v> Jo i- - y dy 

where the neutrino interaction length is defined by Amt = l/{Nji(T^^), the fractional 
energy loss by Ey — E^/ (1 — y) and ^v{E^, 0) — where — ^v{Ev) is related to the 
initial unmodified fluxes as shown in Fig. 1. It turns out to be convenient to solve this 
equation iteratively via the ansatz [35] 

$,(E„X)/$,(E,) = exp [-NAal^,{E,)X] ^.{E,,X) , (2.10) 



d?l!,{E,,X) _ f dy MEy) 



I.e., 



dX l-vME.) 

exp {-Na K(i^.)- <o^(i^.)]X} ^^^^^1^ *.(£;„ X). (2.11) 

(It should be noticed that this solution can be found even more efficiently by using [45] , 
instead of the r.h.s. of (2.10), ex.p[—X/A,,{Ei,,X)] with an effective absorption length 
A,y{Eiy,X) = Aint/ [1 — Zy{Ei,,X)].) According to (2.10), the total shadow factor is given 

by 

S{E,) = — dcp sin ^ exp [-7V^ <j^(£;,)X(^)] *,(£;„ X{9)) (2.12) 

Jo Jo 



where the d(f integration is trivial for our isotropic (anti) neutrino fluxes. The amount 
of material encountered by an upward-going neutrino in its passage through the Earth 
is given by the column depth X{6) which depends upon the nadir angle 6 between the 
normal to the Earth's surface (passing through the detector) and the direction of the 
neutrino beam incident on the detector {9 — 0° corresponds to a beam traversing the 
diameter of the Earth). It is obtained from integrating the density p(r) of Earth [49] 
along the neutrino beam path at a given 9 and is given in Fig. 15 of [20] where X{9) has 
been denoted by z{9). For deflniteness all above formulae have been given for an incoming 
neutrino beam, but similar expressions hold of course for antineutrinos. 

Before evaluating total event rates, we first compare the effects of attenuation and 
regeneration in (2.10) for the initial differential i> + p fluxes in Fig. 1 when they reach 
the Earth's surface {^i,^p{E^) — $(£^,^)/4), with the pure attenuation (absorption) where 
regeneration is omitted = 1 in (2.10)). Some typical representative results for the 
differential shadow factor in (2.10) are shown in Fig. 3 where, for comparison, the pure 
attenuation (absorption) is shown by the dashed curves. The results for u and P fluxes 
are shown separately in Fig. 3(a) and Fig. 3(b), respectively, which are combined in Fig. 
3(c) for the total fluxes in Fig. 1. For the rather flat (initial) cosmic fluxes in Fig. 1, the 
regeneration effect is signiflcant and increases the fluxes at the detector by as much as 30% 
for 9 — 0° and about 10% for 9 — 40° as compared to the purely absorbed (attenuated) 
fluxes. In fact it can even enhance the initial flux at lower energies {E^, ^10^ GeV) at 9 
<, 40° as shown for the AGN-SS flux in Fig. 3. A similar result was originally found in 
[44] investigated in more detail in [35] for (even somewhat flatter) actual neutrino fluxes. 
Below 10^ GeV the shadow factors for antineutrinos in Fig. 3(b) are obviously larger 
than the ones for neutrinos in Fig. 3(a), since a"^ < a^^ , except for the regenerated 
AGN-SS flux at ^ = 0° below 10^ GeV where, in addition, AGN-SS strongly decreases 
for decreasing as can be seen in Fig. 1. At larger energies {Ey > 10^ GeV), where 
the CC and NC cross sections for neutrinos and antineutrinos practically coincide, the 
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shadow factors in Figs. 3(a) and 3(b) become indistinguishable. Clearly, once the energy 

{-) 

is sufficiently high the attenuation factor exp^—NAcr^^^ X) leads to total shadowing as 
seen in Fig. 3 due to the increase of (anti) neutrino cross sections with energy. The smaller 
the nadir angle the lower the energy of complete attenuation. 

For experimental purposes the relevant and interesting quantity is the total u + u flux 
^,^+q{E,^, X) reaching the detector, after its passage through Earth, at different nadir 
angles 9. The four plots in Figs. 4(a) and 4(b) show first the initial u + P fiuxes reaching 
the Earth's surface, ^^+i^{Ei,,X = 0) = $,,+p(£',,) = ^{Ej,)/2 with the original cosmic 
fiuxes $ being given in Fig. 1, and then the fiuxes at the detector X) for the 

three nadir angles 9 = 80°, 40° and 0°. (Recall that 0° is the passage of (anti) neutrinos 
through the center of the Earth.) Due to the large ^u^N cross sections the attenuation 
effects (absorption plus regeneration) reduce dramatically the fiuxes of ultrahigh energy 
(anti) neutrinos, particularly at small nadir angles in Fig. 4. For convenience of reference 
the atmospheric (ATM) neutrino fiux is shown in both Figs. 4(a) and 4(b). The AGN 
fiuxes in Fig. 4(a) stand out above this background ATM neutrino spectrum for neutrino 
energies above about 10^ GeV, whereas the flatter TD and Z-burst fluxes in Fig. 4(b) 
are overwhelmed by the ATM background below 10^ GeV, in particular for smaller nadir 
angles. Nevertheless there remains a window for the observation of cosmic neutrinos (with 

^ 10^ GeV) by underground detection of the energetic when calculating total event 
rates. 

The total upward + /j," event rates are calculated according to (2.1) for the (forth- 
coming) AMANDA-II [2, 50], IceCube [51] and ANTARES [39] detectors with their re- 
spective energy dependent areas A{En) shown in Fig. 5. The required total and flux 
in (2.1) reaching the Earth's surface is given by $i/^+p^ = ^i/^ + = $/2, with the orig- 
inal cosmic neutrino fiuxes $ given in Fig. 1, since the recent discovery of near-maximal 
i/e — and Vn — Vr mixing [52] implies that the originally produced cosmic neutrino 
fiuxes, having v^, : : ratio of 1 : 2 : at some astrophysical source, inevitably 
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oscillate [53, 54] to a ratio of 1 : 1 : 1. The resulting annual total nadir-angle-integrated 
rates are given in Table 1 where the numbers in parentheses are the event rates with no 
regeneration = 1 in (2.10)), i.e. the attenuation of the Earth's penetrating neutrinos 
is just due to naive absorption. Regeneration effects increase the event rates on the av- 
erage by about 20% as compared to the ones where the attenuation of neutrinos is just 
caused by absorption [20, 21, 23]. The contribution to the event rates in Table 1 from 
energies above 10^ GeV becomes, however, negligible and unmeasurably small due to the 
reduction of the initial neutrino fluxes by attenuation with or without regeneration. The 
highest signal rates arise in the AGN models which might be testable for neutrino flux 
energies as large as 10'' - 10^ GeV, i.e. = 10'' GeV. Beyond neutrino energies of 10^ 
GeV, however, present models of cosmic neutrino fluxes are not testable by upward-going 
/i'^ + fj,~ events. Notice that the atmospheric neutrino background becomes marginal 
for neutrino energies above 10^ GeV [20, 21, 35], i.e. Ef"" = 10^ GeV in Table 1, or in 
other words the ATM rate comes entirely from E^, < 10^ GeV. Furthermore the energy 
dependence of the upward-going muon rate in Table 1 will be an important discriminant 
for separating atmospheric and extraterrestrial sources. The nadir angular dependence 
of the annual upward /i'^ + ii~ event rates shown in Fig. 6 for two representative cosmic 
neutrino flux models illustrates in more detail the enhancement of event rates caused by 
regeneration. Apart from the absolute normalizations, these angular distributions are not 
too sensitive to the chosen value of i?™" and clearly favor shallow nadir angles, i.e. large 
9, where the largest amount of events reside. The dashed histograms in Fig. 6 refer to 
events where the neutrino attenuation is caused just by absorption = 1 in (2.12)) with 
regeneration effects disregarded which give rise to the total event rates in parentheses 
in Table 1. Due to maximal i^^ — Vr mixing, the Vt + 9^- flux arriving at the Earth's 
surface may enhance these upward + fi~ rates due to their interaction in Earth via 
Ut-N rX ixX' [55, 56, 57]. Here we disregard such additional contributions since 
Earth-skimming tau-neutrinos [42, 43, 25] will dominate over these interactions and allow 
to test cosmic neutrino flux models at highest energies to which we shall turn later. 
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The upward event rates studied thus far depend obviously rather httle [20] on the spe- 
cific choice of parton distributions since the combination P^S" in (2.3) is rather insensitive 
with respect to different choices: for increasing energies the effect of larger cross sections 
is to increase the probability that a neutrino produces an observable muon, but also to 
increase the attenuation of neutrinos via S on their way to the detector. This is of course 
in contrast to the downward event rates for *'i/^'s that enter the detector from above or 
quasi-horizontally, i.e. where [20] 5" ~ 1. For very high muon-energy thresholds it will 
thus be necessary to observe downward-going muons produced by interactions within the 
instrumented volume. Our expected annual downward /i"*" + /i^ event rates are shown 
in Table 2 which are of course consistent with the ones observed in [21] for appropriate 
fluxes and parameters used there. These results are encouraging and allow to test cosmic 
neutrino fluxes at higher neutrino energies (by about a factor of 10 higher than for the up- 
ward rates in Table 1), provided that downward-going contained events can be observed 
efficiently. For comparison we show again in Table 2 the background atmospheric ATM 
rates which, for the sizeable and measurable cosmic rates, play no role for Ei, 10^ GeV. 
Despite the fact that underground muon detectors (will) have a poor energy resolution, a 
measurement of the energy dependence of the downward-going muon rate will be again 
an important discriminant for separating atmospheric and extraterrestrial sources below 
10^ GeV. (Notice that the downward muon rates in Tables 2 and 3 remain essentially 
unchanged if the detector depth is taken into account. Since underground detectors are 
deployed at a depth of 2 to 3 km, the limited amount of matter above the detector does 
not reduce the shadow factor S{E,y) ~ 1 in (2.3) which holds for all relevant energies in 
Tables 2 and 3. Only at energies as high as Ei, — 10^° GeV it decreases to 5" = 0.99 .) 

Here the question arises whether some other reasonably founded extrapolation to the 
ultrasmall-x region could sizeably enhance our benchmark rates (e.g. in Table 2) calcu- 
lated from cross sections as derived from our nominal QCD-dynamical GRV98 parton dis- 
tributions (or equivalently from the appropriately extrapolated CTEQ3-DIS parametriza- 
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tion). In particular for 'new physics' searches [1, 3, 4, 5, 13, 58, 59, 60] it will be important 
to know to what extent the conventional standard model would allow for (sizeably) larger 
rates. For this purpose wc employ a Regge model inspired small-x description of the DIS 
ep structure function F2^{x, Q^) recently suggested by Donnachie and Landshoff (DL) 
[61]. According to DL, may be written as a sum of three factorized terms fi{Q'^)x~^^ 
to be fitted to DIS HERA data as well as to total photoproduction (Q^ = 0) cross sec- 
tions. This result may be used as a (possibly extreme) guideline for an extrapolation into 
the ultrasmall-x region, dominated by the 'hard pomeron' component Eq ~ 0.4, and in- 
terpolating [62] it smoothly to measured (anti)neutrino structure functions -^2,3 '^(^j Q^) 
by utilizing the CTEQ5 parametrization gives rise to a quicker power-like increase of 
'^vN cross sections at extremely high energies, E^, ^ 10^ GeV, due to the power-hke 
small-x behavior imphed by the DL-fit. This approach will be called [62] DL-I-CTEQ5 
hereafter. The resulting (E^) is shown in Fig. 7 and compared with our nominal 
(radiative GRV98 or CTEQ3-DIS) cross sections as well as with the somewhat smaller 
one corresponding to CTEQ4-DIS. The Regge 'hard pomeron' pole small-x extrapolated 
DL-I-CTEQ5 cross section is, at highest energies, a factor of about 2 larger than our 
nominal cross section predicted by purely radiative QCD-RG evolutions which may be 
considered as a reasonable upper bound implied by a 'conventional' standard model ap- 
proach. Similar remarks hold for the resulting downward fi'^ + /i~ event rates in Table 3 
which should be compared with the nominal ones in Table 2: typically the Regge-inspired 
rates are enhanced by about 20%. It appears that yet much larger rates and cross sections 
than those in Table 3 and Fig. 7 are unlikely to be accommodated by currently conceiv- 
able standard model approaches and, if confirmed by future measurements, might require 
'new' physics ideas for their explanation. 
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3 Event Rates for Surface Telescopes 



For exploring the entire spectrum of cosmic neutrino fluxes up to highest neutrino ener- 
gies of about 10^^ GeV in Fig. 1, large-area ground arrays and/or surface air fluorescence 
telescopes (hke the Pierre Auger Observatory [28, 63] and the Telescope Array [26, 64]), 
where the interaction medium is not the Earth but the atmosphere, will become instru- 
mental since upward-going neutrinos are essentially blocked by Earth for energies above 
10^ GeV. Here, near-horizontal incoming neutrinos will produce electromagnetic and/or 
hadronic extensive air showers which will be observed by Cerenkov radiation and/or fluo- 
rescence detectors. For a, {u + P)-flux reaching Earth, the event rate per second for 
deeply penetrating horizontal showers is given by [28, 26] 

rate [E,, > E.^] = Na Pair / dE.^ dy ^5!d±i^ da(-+-)N(E., y) 

JEth Jo dE^ dy 

where Pair — 10~^ g/cm^ and — for I — e, ji, t and with the original -\- 

flux $ at its production site given in Fig. 1 and the electron neutrino flux has been 
approximated by naive channel counting in pion production and decay {ve/^n = 1/2) 
together with maximal mixing. The relation between the shower energy £'sh and the 
primary neutrino energy depends on the neutrino interaction being considered and 
A{Esh) is the (geometric) detector acceptance. If a detector just measures the shower 
energy (i.e. cannot distinguish between hadronic and electromagnetic showers) we have 
Esh. — Ej, for {ve-\-Ve)N CC interactions, whereas Esh — yE^ for {f^i+V^N CC interactions 
where electromagnetic showers are negligible. This latter relation £^sh = yEy holds also 
for all NC reactions which always produce hadronic showers. If a detector can distinguish 
between hadronic and electromagnetic showers this latter relation remains the same for all 
cases discussed except for {ve-\-Vg)N CC interactions where the neutrino energy has now to 
be shared between the hadronic and electromagnetic shower energy, i.e. E'sh.h = yE^ and 
-£■811,61111 = (1 ~ y)Ey. In Fig. 8 we display the relevant acceptances for the Auger detector 
[28] used for our calculations. We shall not employ the Auger acceptance as calculated 
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by Billoir [65] which has only been presented up to £^sh = 10^*^ GeV since we refrain 
from arbitrarily extrapolating [21] it to 10^^ GeV. The expected annual event rates are 
displayed in Table 4, employing the geometrical acceptances of Fig. 8. Utilizing instead 
the Monte Carlo simulated acceptance in Fig. 8 the resulting total (CC + NC) Auger rates 
would be about a factor of 3 smaller than the ones shown in Table 4. The largest event 
rates arise from (i/g+Pg)-^ CC interactions, for which E'gh ~ E^,. Therefore a measurement 
of the shower-energy dependence will shed some light on the theoretical cosmic neutrino 
flux models in Fig. 1 at ultrahigh energies Ei, > 10^ GeV despite their decrease with 
increasing E^y which, in some cases (TD, Z-burst), is overcompensated by the increasing 
acceptances in Fig. 8. For example, quasi-horizontal air showers with Egh > 10^ GeV will 
test the AGN-M95 flux dominantly at E^ ~ 10^ - 10^° GeV, the TD-SLBY flux at about 
10^ - 10^^ GeV and the Z-burst initiated flux at 10^° - 10^^ Qgy according to the rates 
for the individual E'gh-bins in Table 4. Unfortunately the total amount of air-shower 
events is modest and in one year a few events may be observed in the Auger detector 
assuming modern estimates of AGN-M95, TD-SLBY and Z-burst neutrino fluxes. Most 
other cosmic neutrino flux models, in particular the recently updated and corrected GRB- 
WB flux [11], yield fractions of one event per year and remain undetectable by horizontal 
air showers. This situation remains practically unchanged even for the optimal Regge- 
inspired structure functions resulting in the largest uN cross section in Fig. 7 (event rates 
shown in parentheses in Table 4) which increase our nominal rates by less than 50%. (Our 
results agree of course with the total air shower rates estimated in [21, 28] when using the 
fluxes, structure functions and acceptances employed there.) Our total Auger results in 
Table 4 are compared, where available, with the expected total rates for the TA detector 
[26] which in some cases may be about a factor of 2 larger than for Auger. 

Furthermore, contributions from r-neutrinos, with similar acceptances above 10^ GeV 
as for //-neutrinos [26], may enhance the total air shower rates by as much as 20 to 
30% [26]. We disregard such contributions here, partly because the appropriate detector 
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acceptances for v-,- induced air showers do not exist in the hterature, but mainly because 
Earth-skimming tau-neutrinos [42, 43, 25] will dominate over atmospheric interactions. 



Tau-neutrinos when skimming the Earth [42, 43, 25] will be far more effective in 
producing, besides electromagnetic showers, specific double shower ('double bang') events. 
When ultrahigh energy cosmic neutrinos [Ey ^ 10^ GeV) penetrate and skim the surface 
of Earth they convert to charged leptons which may leave the Earth's surface essentially 
horizontally. For this to happen the neutrinos have to enter Earth near-horizontally at 
some critical nadir angle 0{E^), typically [25, 26] 9{Ey) >, 85", where they travel along 
chords with length of the order of their charged current interaction length A[^f (i?j^). For 
larger nadir angles they rarely interact to produce charged leptons, whereas for smaller 
ones the neutrinos are shadowed by Earth. Electrons lose their energy too quickly in the 
Earth for being detected by ground-level surface telescopes. Furthermore, since ~ 
0.8 X 10~^ cm^/g ^ ~ 6.0 x 10~^ cm^/g in the energy-loss relation (2.6), muons and 
taus can travel 1.5 km and 11 km, respectively, before losing a decade in energy [25, 48]. 
Thus, tau-neutrinos contribute dominantly to Earth-skimming events observed above 
the Earth's surface which consist mainly of electromagnetic showers and of 'double bang' 
showers [25, 26, 42, 43]. The detected tau leptons determine the relevant total event rate 
which is given by [25] 



for an isotropic neutrino flux = = ^/4, according to the definition below Eq. 
(3.1), arriving with nadir angles 9 < 7r/2 and the additional cos^ derives from projecting 
its trajectory (which nearly coincides with the r-trajectory for UHE i/^-'s) onto the surface 
area [66]. The probability (per Er) that a i^t entering Earth with energy Ei, and nadir 
angle 9 produces a r that exits Earth with energy Ej- is 



N^ = 2 dEy dEr d cos 9 K{Ey, 9; E^) 



cos 9 d^y^jEy) 
27r dEy 




(E,)TD 



(3.2) 




2R^ cos 9 



K{E,,9; Er) 



(3.3) 
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with p{r) being the Earth density at distance r from its center (taken from Fig. 14 of 
[20]) with r'^{6,z) = R^^ + — 2R^zcos9, i?© ~ 6371 km, where z is the distance the 
incoming z/^ travels through Earth before converting to a r near the Earth's surface where 
the density is ps — 2.65 g/cm^. Furthermore, CTr = 87.11 /xm and mr = 1777.03 MeV. 
To account for the requirement of clear moonless nights for fluorescence detection, a duty 
cycle of 10% will be assumed in (3.2) as usual, D — 0.1, and the time T an experiment 
runs is taken, for convenience, to be 1 year although, once approved, the experiment 
will run for several years. For reasons discussed above, the integration over the nadir 
angle will be restricted to [25, 26] 90° — ^ < 5° although the results do not depend too 
strongly on this specific choice. The factor of 2 in front of the integral in (3.2) accounts 
for the contribution of incoming antineutrinos since for Ei, ^ 10^ GeV the CC u-,- and 
Pt- cross sections are virtually identical. The effective (geometrical) aperture A^g{ET-) in 
(3.2) as estimated [25] for the TA detector (one station) is shown in Fig. 9 where for 
comparison the one for the far smaller HiRes detector is shown as well. The sensitivity 
of the fluorescence detectors of the Auger Observatory is expected [25] to lie somewhere 
between that of HiRes and TA in Fig. 9, but such a quantity has unfortunately not yet 
been separately published. 

The nominal event rate (based on our GRV98 or CTEQ3-DIS parton densities) for the 
six neutrino sources given in Fig. 1, binned by the r-energy, are given in Table 5 for the 
TA detector assuming 10 observational stations [26]. (We refrain from recalculating the 
rates for the Fly's Eye and HiRes detectors which have been shown to be marginal [25].) 
The expected rates for the largest Regge-inspired uN cross section (cf. Fig. 7) are shown 
in parentheses which are in most cases smaller than the nominal ones (based on GRV98 or 
CTEQ3-DIS in Fig. 7) since the first (attenuation) exponential in (3.3) becomes dominant 
for El, ^ 10^ GeV even for the rather small column depth involved. This implies that 
if event rates are calculated utihzing even smaller UHE uN cross sections than the ones 
in Fig. 7 (due to flatter and thus less steep extrapolations to the ultrasmall Bjorken- 
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X region x ~ M^/s of conventionally fitted parton densities), one easily over estimates 
Earth-skimming rates at highest energies. For comparison the total TA event rates as 
estimated in [26] are, where available, shown as well which are somewhat smaller than 
our results and the ones obtained in [25]. This is presumably due to the fact that the 
realistic TA aperture is expected to be somewhat smaller [67] than the one calculated in 
[25] and shown in Fig. 9 because of effects such as the detector response, light propagation 
in air and night sky background. The total rates for the Auger detector are also shown 
in Table 5 which we have calculated for our neutrino fluxes in Fig. 1 using the effective 
Auger aperture [43] where, besides the dominant electromagnetic showers, multi-bang 
events have been also taken into account for the Monte Carlo simulations (for definiteness 
we used the average aperture in Fig. 9 of [43] as obtained from the 'BS+PP+DIS low' 
continuous energy loss model) . In any case, in contrast to the quasi-horizontal air showers 
in Table 4, the shower event rates produced by Earth-skimming i/^-'s are sufficiently large 
for testing and exploring the large energy tail 10^ GeV) up to 10^^ GeV of most cosmic 
neutrino sources in Fig. 1, except for the recently updated and corrected GRB-WB flux 
[11] for which the rates remain prohibitively small. 

In order to learn to what extent the shower rates measured in certain bins of Et- as 
shown in Table 5 can explore and test the energy profile of the theoretical cosmic neutrino 
fluxes ^{E^) in Fig. 1, we display in Table 6 the individual £'i,-binned contributions to 
our nominal TA-rates for the £^T--bins given in Table 5 for three representative neutrino 
fluxes. Apart from the less steep falling Z-burst flux in Fig. 1 which contributes over 
a wider E^-range for a given £^T--bin, measured showers with E^ = 10® — 10^ GeV will 
test cosmic neutrino fluxes mainly within E,y — 10® — 10^° GeV, whereas increasing the 
shower energy to E^- — 10^ — 10^° GeV, for example, will dehneate a similar £',^-bin. 
El, — 10^ — lO^'^ GeV. The same holds true for even higher values of Et- Therefore from 
an experimental point of view it will be important to measure the detailed energy (E^-) 
dependence predicted in Table 5 for exploring cosmic neutrino sources at highest energies. 



17 



4 Implications of Low Scale 'Extra Dimensions' 



Having calculated so far in the previous Section the muon and shower event rates caused by 
extragalactic neutrinos according to conventional (nominal) and optimal (Regge inspired) 
standard model expectations, we will now finally give for comparison a few examples 
of distinctively different event rates which derive from 'new' physics ideas. Although 
highly speculative and by far non-unique, let us choose string theories with large 'extra 
dimensions' 6 where only gravity propagates in the 4 + 5 dimensional bulk of spacetime 
(a recent brief review can be found in [68]). To obtain the effective low-energy theory in 
3+1 dimensions, these extra dimensions are compactified to a common scale R which 
may be relatively large, corresponding to a small scale 1/i? of new physics, and is related 
to standard Newton gravity (Planck scale) via Gj^^ ~ Mp'^R^ with Ms ~ TeV being the 
low energy effective string scale [69]. (For example, R is of the order of 1 mm for 6 = 2.) In 
addition to the usual Standard Model (SM) particles, an infinite tower of massive Kaluza- 
Klein (KK) excitations appears in the effective four-dimensional theory, corresponding 
to the massless graviton in 4 + 5 dimensions, and Ms plays the role of an ultraviolet 
cutoff for the summation over all relevant spin-2 KK states. Such a scenario opens up 
the interesting possibility of massive spin-2 t-channel exchange which results in a more 
rapid growth with energy as compared to SM cross sections derived from massive vector 
boson exchange considered thus far. Being flavor-neutral this 'graviton' KK exchange {G) 
contributes only via NC processes to the total SM uN cross section. The cross sections for 
the partonic subprocesses i/V— > i^V and vg ^ vg bs mediated by G read [58, 59, 70, 69] 

^ = ^|^[325^ + 64J^t + 42^^t^ + 105f^ + tV(.;t) (4.1) 

^ = ^^[2^' + 4J'^^+3^^P + ^f^]F(J,£) (4.2) 

s 

ior 5 — 2 extra dimensions and the terms in square brackets are symmetric under s <-> it. 
The function F{s,i) refers to some unitarization (interpolation) procedure [58, 59, 71] in 
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order to extrapolate the perturbative predictions at s ^ M| in (4.1) and (4.2) to s ^ M|. 
We shall use two representative choices [58], 

F(s,t) = ^ ri + 101n(l + s7M|)l (4.3) 

^ ^ (M| + s)2(M|-i) L \ / sJi \ J 

(M2+s)(M|-t) 

as well as choosing Ms — 1 TeV for our subsequent quantitative analysis. The contribution 
to the total SM uN cross section due to 'graviton' KK exchange is thus given by 

f=<l,<l,9 

with s — xs, t — t = —Q^ — —xys and only the dominant light quarks q — u, d, s are 
taken into account. The same expression holds for PN scattering. In Fig. 10 our nominal 
SM total (CC + NC) vN cross section is compared with the additional KK contribution 
where the dominant gluon initiated component in (4.5) is, at highest energies, about a 
factor of 2 to 4 larger than the sea contributions of (anti)quarks depending on whether 
the interpolation (4.4) or (4.3), respectively, is used. (It should be noticed that these 
cross sections are in agreement with the upper bound derived in [72]. A more recent 
updated estimate [73], however, being based on the full exposure of the AGASA and 
Fly's Eye experiments as well as on a larger cosmogenic neutrino flux [74] together with a 
larger nucleon source cutoff energy, would rule out the largest strongly rising cross section 
(dashed curve) in Fig. 10 at highest neutrino energies.) These total SM + KK uN cross 
sections can rise to hadronic mb-scale values [58, 59] at highest neutrino energies and 
the resulting annual total event rates for the relevant quasi-horizontal air showers, shown 
in Table 7, are enormous when compared with the small SM rates in Table 4. Even for 
highest neutrino energies (i.e. highest threshold energies in (3.1)) the number of events 
remains large, except for the steeply falling AGN-SS and GRB-WB fluxes in Fig. 1. It 
should be easy for the future Pierre Auger Observatory to observe such dramatically large 
hadronic shower event rates as implied by low-scale large 'extra dimension' scenarios. 
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(Given the uncertainties in the cosmic neutrino flux predictions, the eventual detection 
of such event rates wiU be subject to interpretation, because they are proportional to the 
product of the flux and the cross section. Ways to discriminate between large neutrino 
cross sections and large neutrino fluxes have been suggested recently by looking at the 
zenith angle distribution of the neutrino induced shower events, i.e. by comparing the 
energy spectrum of upgoing and nearly horizontal to downgoing shower events [66, 75]. 
Such a measurement may at least in principle allow to determine a^^^^^^ independently 
of the neutrino flux.) The viability of such 'new' physics scenarios may be even tested 
by the smaller ground level AGASA [76] and Fly's Eye (HiRes) [77] observatories: using 
their estimated shower acceptances [78, 77], the expected rates for AGASA are about a 
factor of 25 smaller than the ones in Table 7, and a similar reduction holds for the Fly's 
Eye (HiRes) observatory. 



5 Summary 



In order to explore and test the sources of cosmic UHE neutrinos, the calculation of the 
relevant '^vN cross sections requires the knowledge of parton distributions /(x, Q^), with 
f — 0.1 Q.-I 9 1 ultrasmall Bjorken x ~ M^^/ below the range constrained by present ex- 
periments, and at scales = M^. Highest neutrino energies Ej, — s/2M]sr ~ 10^^ GeV 
require extrapolations down to x ~ 10~^. Within standard QCD RG-evolutions, such 
extensive small-x extrapolations can be uniquely and reliably predicted by the dynam- 
ical (radiative) parton model [32] which proved to provide reliable high energy small-x 
predictions in the past. (Interestingly, consistent BFKL model resummations in leading 
order QCD yield remarkably similar results [35].) All resulting relevant ^u^N cross sections 
turn out to have a typical uncertainty of about ± 20% at highest neutrino energies of 
10^^ GeV [36]. We have adopted this 'nominal' approach for all our calculations which 
are based on the canonical radiative GRV98 parton distributions [32] . Similar results are 
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accidentally obtained from the conventionally fitted CTEQ3-DIS parametrizations [37] 
at x ^ 10~^ with their assumed fixed-power extrapolation to x < 10~^, whereas the 
CTEQ4-DIS parametrizations [38] somewhat underestimate the neutrino-nucleon cross 
sections by about 20% at Ei, ~ 10^° — 10^^ GeV. In calculating event rates caused by cos- 
mic UHE neutrinos we have also employed recent Regge-model inspired small-x structure 
functions [61, 62] in order to learn about possible maximal rates which may be accommo- 
dated by 'conventional' standard model approaches and which do not necessarily require 
'new' physics interpretations. 

For the representative set of modern cosmic neutrino fluxes in Fig. 1 we flrst analyzed 
/i"*" + n~ event rates caused by upward-going UHE neutrinos for modern (future) under- 
ground detectors (such as ANTARES, AMANDA II and IceCube), with special emphasis 
on NC regeneration effects which populate the lower energy part of the attenuated flux 
spectra. Regeneration increases the non-regenerated event rates [20, 21, 23] on the av- 
erage by as much as 20%. In particular AGN models, yielding the highest signal rates, 
might be testable in this way for neutrino energies as large as 10^ — 10^ GeV, i.e. well above 
10^ GeV where the atmospheric neutrino background becomes marginal. In addition we 
demonstrated in detail in Table 1 to what extent high-E^;^ bins of the original fluxes in 
Fig. 1 remain experimentally accessible after their depletion and regeneration due to NC 
interactions. Measurements of these energy-dependencies will be crucial for exploring the 
nature of various cosmic neutrino sources. 

In contrast to these upward /i'^ + ii~ rates, the downward-going muon rates depend 
strongly on the speciflc choice of parton distributions and would allow to test cosmic neu- 
trino fluxes at higher E^y (about a factor of 10 higher than for the upward rates, cf. Table 
2), provided that downward-going muons, produced by interactions within the instru- 
mented volume, can be observed efficiently. The largest Regge-inspired cross section (s) 
in Fig. 7 give rise to downward -|- event rates (Table 3) which are about 20 — 30% 
larger than our nominal ones in Table 2, and may be considered as reasonable upper 
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bounds of present standard model approaches. 

To explore the entire spectrum of cosmic neutrino fluxes up to highest neutrino energies 
of about 10^^ GeV, one has to resort to quasi-horizontal air showers and Earth-skimming 
i^^'s to be detected by large area ground arrays and surface air fluorescence telescopes (such 
as the Pierre Auger Observatory and the Telescope Array), since upward-going neutrinos 
are essentially blocked by Earth for energies above 10^ GeV. The detailed shower energy 
Esh dependence, which practically coincides with the one of the incoming cosmic neutrinos, 
of the annual event rates for the Auger detector has been presented in Table 4 as well as 
the results derived from the largest Regge inspired cross sections which can be at most 50% 
larger than our nominal rates. The most efficient way, however, to test cosmic neutrino 
fluxes at highest energies, 10® ^ Ey ^ 10^^ GeV, will be Earth-skimming tau-neutrinos 
[42, 43, 25] which, when converting to r-leptons that escape Earth, produce dominantly 
near-horizontal electromagnetic showers as well as specific 'double-bang' events. Our 
nominal rates are presented, for various bins in in Table 5 and compared with the 
expected rates for the Auger detector. Here the largest Regge-inspired vN cross section 
(cf. Fig. 7) implies in most cases smaller rates than our nominal ones (based on GRV98 or 
CTEQ3-DIS) since the attenuation (absorption) exponential in (3.3) becomes dominant 
at neutrino energies above 10^ GeV despite the rather small column depth involved. In 
order to learn to what extent the measured shower rates in certain bins of E^. will be able 
to delineate specific features of the energy profile of the various theoretical cosmic neutrino 
fluxes ^{Ey) in Fig. 1, we presented in Table 6 the individual i^j^-binned contributions to 
the nominal TA~rates for given i?T--bins. In any case it will be important to measure the 
detailed energy {Egh and Et) dependencies for exploring extragalactic neutrino sources at 
highest energies of 10^^ eV. 

Finally, as an illustration of possible implications of 'new' physics, we calculated highly 
speculative expectations of string scenarios with large 'extra dimensions' at low TeV- 
scales. The additional contributions, due to the exchange of massive 'graviton' spin-2 
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Kaluza-Klein excitations in neutral current neutrino-quark (gluon) scattering to the rele- 
vant quasi-horizontal hadronic air showers may be enormous as compared with the mod- 
erate SM event rates for the Auger detector and could be detected in the not too distant 
future. The viability of such 'new' physics scenarios may be tested even by the smaller 
ground level AGASA and Fly's Eye (HiRes) detectors. 
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Table 1: Total upward fj,~^+fi~ event rates per year from j/^N and f'^N interactions in rock calculated according 
to (2.1) for various muon energy thresholds E™*" and the appropriate cosmic neutrino fluxes in Fig.l. The 
numbers in parentheses refer to the event rates with no regeneration, i.e. ^ =1 in (2.10). The background ATM 
rates are shown for comparison. 
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Tabic 2: Total downward /i+ + event rates per year arising from Vf^N and Vf^N 
interactions in a detector with an effective volume 14// = ^e// x 1km = Ikm^ of 
water calculated according to (2.3) with 5 = 1 for various muon energy thresholds 
E™^. The background ATM rates are shown for comparison. 
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Table 3: Total downward /i+ + n~ event rates per year for a detector with 
Veff = Ikm'"* as in Table 2 but using the DL+CTEQ5 structure functions [62] 
with their Regge model inspired small-x extrapolation [61]. The background ATM 
rates are shown for comparison. 
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Table 4: Total annual event rates for the Pierre Auger Observatory for horizontal air showers induced by {ve + f'e)]^ 
and + S'y:t)N CC and NC interactions calculated according to (3.1) using the geometrical acceptances of Fig. 8 and 
the appropriate cosmic neutrino fluxes of Fig.l. The numbers in parentheses refer to the rates as obtained by using 
the (optimal) DL+CTEQ5 structure functions [62] with their Reggc model inspired small-x extrapolation [61]. For 
comparison the expected rates for the Telescope Array are also displayed, where available, which have been estimated 
[26] using CTEQ4-DIS parton distributions and assuming a duty cycle of 0.1 for observing the atmospheric fluorescence 
light. 
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Table 5: Annual nominal event rates produced by Earth-skimming z/^'s for the 
fluorescence Telescope Array detector (10 stations) calculated according to (3.2) 
with D=0.1 and the aperture for the detections of r leptons through their dom- 
inant decay to electromagnetic showers taken from [25]. The numbers in paren- 
theses are the rates resulting from the largest Regge-inspired t'N cross section in 
Fig. 7. The estimated total TA rates of [26] arc also displayed, where available. 
The total Auger rates are obtained from folding the effective Auger aperture [43] 
('BS+PP+DIS low' curve in Fig. 9 of [43]) with the appropriate cosmic z/T-fluxes 
in Fig.l. These latter event rates refer to the Auger ground array (not the fluo- 
rescence detector) with a duty cycle of 1. 







detectors 




flux 


E (GeV) 


TA 


TA([26l) 


Auger([43]) 




1 ns _ 1 n9 


1 07 


fl 19'l 








J-U J-\J 


OOOfi 


(0 00061 








in^o - mil 

IQii - 10^2 




[-) 

(-) 










1.07 


(1.12) 




1.15 




10^ - 10^ 


47 4 


(Af> 1 








1 - 1 n^o 

lU — lU 




fl 9 R\ 
±z,.o I 








inio _ inii 

1 U -LU 

IQii - 10^2 


0.25 


(0.19) 

K J 








tot^i 1 


62.5 


(59.1) 


42.8 


72.8 




1 — 1 


n nns^ 

\J .yjyjoo 


fn no36'i 

U .UUOU I 








-LVJ -LVJ 


0002 

yj .yjyjyj ^ 


fn 0002^ 






OR R-WR 


inio - inii 
10" - 10^2 




(-) 

(-) 








total 


0.0037 


(0.0038) 




0.0039 




10^ - 10^ 


7.62 


(7.22) 








109 - lO^o 


5.21 


(4.35) 






TD-SLRY 


10^° - 10" 
10" - 10^2 


0.54 

0.008 


(0.40) 

(0.005) 








total 


13.4 ( 


11.98) 




17.0 




10^ - 109 


0.19 


(0.17) 








109 - 10^° 


0.23 


(0.18) 






TD-SLSC 


lO^o - 10" 
10" - 10^2 


0.05 
0.002 


(0.04) 
(0.001) 








total 


0.47 


(0.39) 


0.34 


0.66 




10^ - 109 


0.34 


(0.29) 








109 - 10^0 


0.80 


(0.63) 






Z-burst 


lO^o - 10" 
10" - 10^2 


0.27 

0.017 


(0.19) 

(0.011) 








total 


1.43 


(1.12) 


1.19 


2.18 



Table 6: Contributions to the nominal TA event rates for given Er-hins 
in Table 5 from various increasing ii^^-bins of three representative cosmic 
neutrino fluxes. 
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Table 7: Total annual event rates for the Pierre Auger Observatory for hori- 
zontal air showers as implied by the SM and the additional contributions due 
to the exchange of KK gravitons in a low-scale scenario with large 'extra di- 
mensions'. The rates are calculated according to (3.1), with the additional 
graviton KK exchange contribution given in (4.5), and using the geometrical 
acceptances of Fig. 8. The first entry corresponds to using the unitarity inter- 
polation (4.3) (dotted curve in Fig. 10) and the second entry is the result of 
using (4.4) (dashed curve in Fig. 10). 
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Figure 1: Representivc differential fluxes of muon neutrinos (y^i+v^) from ac- 
tive galactic nuclei (AGN-SS [9] and AGN-M95 [8]), gamma ray bursts (GRB- 
WB [11]), topological defects (TD-SLSC [14] and TD-SLBY [15]) and Z-bursts 
[18]. Due to naive channel counting in pion production and decay at the pro- 
duction site [ve : fix = 1 : 2) and maximal mixing, : : z^r = 1 : 1 : 1, 
these fluxes are divided equally between ^, e and r neutrinos when they reach 
Earth (i.e. will be devided by a factor of 2). The background angle-averaged 
atmospheric (ATM) neutrino + flux [19,20,21] is shown for illustration 
by the long-dashed curve. 
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Figure 2: Mean ranges in rock of muons produced in CC interactions of neutri- 
nos with energies according to (2.5) with E™" =lTeV. The analytic ranges 
follow from (2.7) for constant values of and and the numerical range fol- 
lows from (2.6) using (2.8). The Lipari-Stanev range ist taken from [48] which 
extends only up to 10^ GeV. 
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Figure 3(a): The differential shadow factor in (2.10) for three different nadir 
angles of cosmic neutrino fluxes incident at the Earth's surface with an initial 
flirx ^„{E,y,X = 0) = ^„{E^) = ^{E^)/4 with the original total + flux 
^{E^) being given in Fig.l. The dashed curves describe the attenuation just 
due to absorption where regeneration is omitted, * s 1 in (2.10). 



I lllllll — I I I llllll — I I I llllll — I 1 1 llllll — I 1 1 llllll — I 1 1 llllll — I 1 1 llllll 

AGN-SS 



TTT1I| 1 I I lllll| 1 I I lllll| 1 I llllll| 1 I I lllll| 1 I I lllll| 1 I I lllll 

AGN-M95 




10^ 10^ 
Ev [GeV] 

Figure 3(b): As in Fig. 3(a) but for antineutrinos where the initial flux at the 
Earth's surface is again given by ^p{Ei,,X = 0) = ^p{Ei,) = $(i?i.)/4 with $ 
being given in Fig.l. 
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Figure 3(c): As in Fig. 3(a) but for the total neutrino plus antineutrino fluxes 
where the initial total flux at the Earth's surface is given by ^^+p{E^,X = 
0) = ^„+p{E„) = ^{E„)/2 with the original $ given in Fig.l. 
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Figure 4(a): The initial Vf_i+v^ flux at the Earth's surface <I>^-|_p(E^,X=0)= 
$i,+p(Ei.)=$(Ei,)/2, with $(Ej.) in Fig.l, and the fltix at the detector <^„^o{^^,^) for 
three different nadir angles corresponding to two models for AGN neutrinos, AGN-SS [9] 
and AGN-M95 [8], and to neutrinos from gamma ray bursts, GRB-WB [11]. The division 
by a factor of 2 of the original cosmic fluxes $(Ejy) in Fig.l is due to maximal mixing. 
The background atmospheric (ATM) neutrino flux is also shown for comparison. 
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Figure 4(b): As in Fig. 4(a), but showing neutrino fluxes from two topological defects 
models, TD-SLBY [15] and TD-SLSC [14], and from Z-bursts [18]. 



LU 
< 



0.1 



















IceCube 




AMANDA-II 




ANTARES 



0.01 

10" 10* 10° 10° 10' 10° 

E^[GeV] 

Figure 5: Effective energy dependent areas for the underground detectors 
ANTARES [39], Amanda-II [50] and IceCube [51] used for our calculations 
of total event rates. 
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Figure 6: Nadir angular dependence of upward /i+ + fjr event rates per year 
for two representative neutrino flux models in Fig.l. The dashed histograms 
refer to events where the neutrino attenuation is caused just by absorption 
(^f = 1 in (2.12)) with no regeneration. For definiteness an effective energy- 
independent detector area of Ae/f = 1 km^ has been used in (2.1) (cf. Fig. 5). 
The TD-SLBY rates have been multiplied by a factor of 200. 
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Figure 7: Total (CC+NC) j^N cross sections. The dasliecl-dotted curve refers 
to the cross section calculated using the DL+CTEQ5 structure functions [62] 
with their Regge model inspired small-x extrapolation [61]. Our nominal re- 
sults are based on GRV98 [32], or equivalently on the CTEQ3-DIS [37] parton 
distributions. For illustration the results based on CTEQ4-DIS [38] are shown 
as well. 
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Figure 8: Acceptances versus shower energy needed for calculating near- 
horizontal air shower event rates according to (3.1) for the Auger detector. 
The solid curves correspond to a geometrical integration for electromagnetic 
and hadronic showers of Capelle et al. [28] and the crosses are the results of a 
Monte Carlo simulation of all showers [28] which can be considered as a lower 
bound of the acceptance. For comparison the acceptance for all showers as 
estimated by Billoir [65] is shown by the dotted curve. 
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Figure 9: Effective (geometric) aperture estimates [25] for the detection of 
Earth-skimming t leptons through their most promising and dominant decays 
to electromagnetic showers for the Telescope Array (one station) and the HiRes 
detector. The TA aperture is exptected to be somewhat smaller [67] due to 
effects such as the detector response, light propagation in air, and night sky 
background. 
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Figure 10: The nominal total vN cross section in the SM(CC+NC) compared 
to a scenario with large 'extra dimensions' using two different unitarity extrap- 
olations between perturbative and nonperturbative regimes: the dottet line 
refers to Eq. (4.3) and the dashed one to Eq. (4.4). 



